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We study the quantum dynamics of a binary mixture of Bose-Einstein condensates (BEC) in a 
double- well potential starting from a two-mode Bose- Hubbard Hamiltonian. Focussing on the regime 
where the number of atoms is very large, a mapping onto a SU{2) spin problem together with a 
C/2 ' Holstein-Primakoff transformation is performed. The quantum evolution of the number difference 

of bosons between the two wells is investigated for different initial conditions, which range from 
the case of a small imbalance between the two wells to a coherent spin state. The results show an 
instability towards a phase-separation above a critical positive value of the interspecies interaction 
C . while the system evolves towards a coherent tunneling regime for negative interspecies interactions. 

A comparison with a semiclassical approach is discussed together with some implications on the 
experimental realization of phase separation with cold atoms. 



PACS numbers: 03.75.Lm, 67.85. Fg, 74.50.-|-r 



I. INTRODUCTION 



Bose-Einstein condensates of dilute, weakly interacting gases offer a unique possibility for exploring many-body 
dynamics, the role of quantum fluctuations and in general macroscopic quantum coherence phenomena thanks 
to a wide tunability of the interaction parameters. Indeed several experimental strategies can be devised in order 
to pursue this task, which range from the direct control via magnetic Feshbach resonance techniques [2| to the 
\^ • transverse confinement in a quasi one dimensional system Q as a way to increase the inter-atomic interaction. 
0^ , Finally, the introduction of an optical lattice whose depth can be tuned allows one to decrease the kinetic term in the 
■"nI" ■ Hamiltonian. Within the tight binding approximation such systems are described by the Bose-Hubbard Hamiltonian, 
] whose parameters are the hopping frequency Ej between nearest neighbor lattice sites, the onsite interaction strength 
. Ec and the total atoms number N. When the ratio exceeds unity, a quantum phase transition from a superfluid 

, to a Mott insulator Q takes place and the system enters a quantum regime characterized by strong correlations. The 
. . ' simplest Hamiltonian of this kind that one can devise is the Bose-Hubbard dimer [Bj, which describes the physics of 
^ ] two weakly coupled condensates. It can be mapped onto a SU{2) spin problem and is deeply related to the physics 
of Josephson junctions Furthermore, if the mean field approximation is considered one obtains the Gross- 

Pitaevskii theory which gives rise to a variety of phenomena, ranging from Josephson oscillations Q to macroscopic 



. quantum self-trapping (MQST) and ac and dc Josephson like effect all experimentally observed in the last 
decade 

More recently, after the experimental realization of two-species BECs [3l [lH [3 , the theoretical analysis on weakly 
coupled condensates has been successfully extended to a binary mixture of BECs in a double well potential [l3| [i3 [Toj 
(20||2l|[2^. The semiclassical regime in which the fluctuations around the mean values are small has been deeply 
investigated and found to be described by two coupled Gross-Pitaevskii equations. By means of a two-mode ap- 
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proximation such equations can be cast in the form of four coupled nonhncar ordinary differential equations for the 
population imbalance and the relative phase of each species. The solution results in a richer tunneling dynamics j23j . 
In particular, two different MQST states with broken symmetry have been found [131, where the two species localize 
in the two different wells giving rise to a phase separation or coexist in the same well respectively. Indeed, upon 
a variation of some parameters or initial conditions, the phase-separated MQST states evolve towards a symmetry- 
restoring phase where the two components swap places between the two wells, so avoiding each other. Furthermore 
the coherent dynamics of a two species BEC in a double well has been analyzed as well focussing on the case where 
the two species are two hypcrfinc states of the same alkali metal p^ . 

In a recent paper [2^ we studied the quantum behaviour of a binary mixture of Bose-Einstein condensates (BEC) 
in a double-well potential starting from a two-mode Bose-Hubbard Hamiltonian. We analyzed in detail the small 
tunneling amplitude regime where number fluctuations are suppressed and a Mott-insulator behaviour is established. 
Within this regime we performed a perturbative calculation up to second order in the tunneling amplitude and found 
the stationary states. In order to carry out analytical calculations we focused on the symmetric case of equal nonlinear 
interaction and equal tunneling amplitude of the two species. Furthermore we restricted to the case in which the two 
species are equally populated and imposed the condition of equal population imbalance of the species a and b between 
the two wells. Then, the dynamics of the junction was investigated in correspondence of a completely localized 
initial state. In order to avoid the above restrictions on the parameters range, here we focus on the two-mode Bose- 
Hubbard Hamiltonian describing the two-species BEC (a and b) in a double well when Ng, Nb ^ 1, and perform a 
mapping onto a SU{2) spin problem together with a Holstein-Primakoff transformation l26l 1271 . As a result we obtain 
a Hamiltonian of two decoupled quantum harmonic oscillators, similar to that of Ref.[28|, whose stationary states 
are readily found. The quantum evolution of the number difference of bosons between the two wells is investigated 
in detail in correspondence of a variety of initial conditions, which range from an initial state with small imbalance 
between the two species to a coherent spin state. The whole parameters space is explored by tuning the population, 
the tunneling amplitude and the nonlinear interaction for each species as well as the interspecies interaction in a 
wide range, from a symmetric to a strongly asymmetric case. Finally a detailed comparison with a semiclassical 
approach is given. Let us notice that Holstein-Primakoff transformation makes the system exactly solvable in the 
weakly interacting regime of interest in this work and that simplifies the study of the tunneling dynamics as well as 
the phase separation phenomenon. This is the main advantage of the approach chosen. 

The paper is organized as follows. In Section 2 we introduce the model Hamiltonian we study within the two mode 
approximation. A Holstein-Primakoff transformation is performed and the semiclassical limit is taken followed by a 
decoupling of the bosonic degrees of freedom for each species. As a result the Hamiltonian can be rephrased in terms 
of two independent harmonic oscillators, whose stationary states are derived in Section 3. In Section 4 the quantum 
dynamics of the system is discussed in correspondence of two different initial conditions: small imbalance between 
the two wells and coherent states. A wide range of values of interspecies interaction is explored and the crossover to 
an unstable regime with phase separation is found. In Section 5 the classical equations of motion are derived and a 
comparison with the quantum counterpart is carried out. Finally some conclusions and perspectives of this work are 
briefly outlined. 



II. THE MODEL 

A binary mixture of Bose-Einstein condensates [l^ [l^l loaded in a double- well potential is described by the Hamil- 
tonian H = Ha + Hb + Hab where: 

H, = I ' 



Hab = 9ab / dl^lp^lp'^lpa'^'b- (2) 



Here gu = °" is the intraspecies coupling constants, being nii the atomic mass and an the s-wave scattering 
lengths; gab = "^^^ is the interspecies coupling constant, where niab — m°+mi, ^^'^ reduced mass; Vi (l^) is the 
double well trapping potential and, in the following, we assume Va (T^) = Vb {l^) = V ('P'); (T^) , i^i (^), i ^ a,b 
are the bosonic creation and annihilation operators for the two species, which satisfy the commutation rules: 

, )] = K (^) . (^')] - 0, (3) 
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and the normalization conditions: 

j dl^\^,{-t)\' = N,- i = a,b, (5) 
Ni, i — a,b being the number of atoms of species a and b respectively. The total number of atoms of the mixture is 

N = Na+Nb. 

A weak link between the two wells produces a small energy splitting between the mean-field ground state and the 
first excited state of the double well potential and that allows to reduce the dimension of the Hilbert space of the 
initial many-body problem. Indeed for low energy excitations and low temperatures it is possible to consider only such 
two states and neglect the contribution from the higher ones, the so called two-mode approximation [2^ QfsO]. In 
this approximation the Hamiltonian ([1]) can be written in terms of the the annihilation operators, = (ig + le), 

qr = {ug — tte) and = {bg + b^), bu ^ ^ i^g ^ ^e) and the corresponding creation operators, where ag,ae 
and bg, be are the annihilation operators of a particle in the ground and in the first excited state. 
When introducing the angular momentum operators: 

4^'^{b+bL + btbi,), j'y^^{b+bL~btbR) j',^'^{b+bn-btbL), 

where the operators Jf , jf , i = x,y, z, obey to the usual angular momentum algebra together with the relation: 

(J-^f = ^(^ + 1), (J'')^ = f (f + 1), (7) 
the Hamiltonian of the double species Bose-Josephson junction can be written in the form: 

+^abJzJz ~ DabJxJx- (8) 

where Ka^b are the tunneling amplitudes between the two wells, A(j,6, Aqj, are the intra- and interspecies interactions 
respectively, while the terms Ca and Da.b describe two-particle processes [2^. The form ([8]) was previously discussed 
in the classical limit in [23] , where it was shown to be lead to equations of motion equivalent to the Gross-Pitaevskii 
equations. For Aab = Dab = in Eq. ([5]), the Hamiltonian reduces to a sum of two Lipkin-Meshkov-Glick (LMG) 
model [111, H^l Hamiltonian, one for each species. For Aab 7^ or Dab 7^ 0, the two LMG models are coupled. Within 
the experimental parameters range it is possible to show that d <s:_ Ai, Ki, i ~ a, b, and Dab ^ Aab 2^ 6][23|, then in 



the following we put Ca = Cb = and Dab = 0, which corresponds to neglecting the spatial overlap integrals between 
the localized modes in the two wells. In this way the binary mixture of BECs within two-mode approximation maps 
to a two Ising spins model in a transverse magnetic field, whose Hamiltonian is: 

H=l^a {J^f - KaJi + i A, {Jl) ' - K,JI + Aabr.Jl (9) 

Let us briefly discuss the symmetries of the Hamiltonian 0. First, when Aab = 0, the Hamiltonian decouples into 
H ~ Ha 4- Hh, and e"'^= H^ie~^'^^^ = H^i for v, v' G {a, 6}. Therefore, the eigenstates of H can be sought in the form 
of eigenstates of e""^^ and e""'»=. Since e^"''^ = e*'^^", these eigenvalues are ±1 when N^, is even, and ±i when N^, is 
odd. So the Hilbert space breaks down into four sectors indexed by the eigenvalues of e'^'^^^ and e"'^^ . Then, turning 
on Aab 7^ 0, the only remaining symmetry is e^^'^''^^'^''^ H^ie~^^'^'^^^'^^\ so that only two independent sectors remain. 
These sectors are formed by the combination in pairs of the four sectors obtained for Aab = 0. 
Let us now make the rotation: 

P — r 

, z^a,b. (10) 



To proceed we perform the Holstein-PrimakofF transformation |26l - l28| in order to map the angular momentum opera- 
tors into bosonic ones and focus on the regime with large number of atoms Na, Nb ^ 1 and weak scattering strengths 

Ka{b) > Aa, Ab, Aab: 

J? = J" - flta = - b^b 



Jl = V2 J° - gtaa , ,/| = V2 ^ b^bb , (11) 
= atV2 J° - ato = bW^Jb - b^b 
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where = J^zL iJ'L, i = a,b, J' = Ni/2 i = a, b, thus leading to the Hamiltonian: 



2 J" (J" + 1) - 2(J° - a^a)^ + \J - aU){2J'' - 1 - a^a)a^ + {a^f^{2J'' - ata)(2J« - 1 - a^a) 



' 8 

Agb 
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2J''{J'' + 1) - 2(J'' - 6^6)2 + ^{2Jb -b'<b){2Jb -l-b'tb)b^ + (6t)2y^(2J« - 6t5)(2J« - 1 - 5t5) 



y^(2J'' - ata)(2J'' - 6t6)a6 + fetflt y^(2J'» - ata)(2J'' - b'<b) 
+ \/2Jb - b^ba^b\/ 2 J« - a^a + \/2J'^ - ab^ a\/2J'' - ata 



(12) 



Here a and b are boson annihilation operators for each species. In this representation, the operators e*'^'^= are equal to 
gnT[j"-u''v) g^j-^j their action is simply v —v. For Kab ~ 0, the Hilbert space of (|12p thus breaks into four different 
sectors, according to the parity of a^a and while for Kab 7^ 0, it breaks into two different sectors depending on 
the parity of a^a + b^b. The physical Hilbert space is restricted to < a^a < Na and < &t& < Nb. 

Since we are considering a large number of atoms, we have J°',J^ 3> 1, while the condition Kg_(b) ^ Aa,Af,,Aa6 
implies (a+a) ^ 2J" and (&"''&) ^ 2J''. Under these assumptions one can use the linearized Holstein-PrimakofF 
transformation [2^ (i.e. J| — J' — s^s, Jf = \/2 J^s, Jf = s^V2 J* with s ~ a,b) and derive the effective Hamiltonian: 



a + a+ \ / a + a+ 



2AabVJ^ 



a+\ fb + b 



KaJ" - Kb.r + /v^a+a + A'66+6. 



(13) 



In order to decouple the degrees of freedom of each bosonic species let us introduce the following harmonic oscillator 
coordinates and momenta, qi, pi, i = a,b: 



Pa^^ia-a+) p,^^{b-b+)' 
which satisfy the usual commutation rules [qi,pj] = i,j ~ a,b. Then, by defining: 



(14) 



— _2a 



_ qb 



Pa = VKaPa, Pb = VKbPb, 

(where [Qi,Pj] = iSij, i,j — a,b) and, by dropping constant terms, Eq. (|13p can be written in a matrix form as [2 



(15) 



H2B.J.J — 2 



Q^cD^Q + P'P 



(16) 



where 



^ab ^b 



(17) 



and = {Qa,Qb), P^ ~ {Pa,Pb) (the symbol stands for the transpose); ijjf ~ KiJ'^Ki + Kf, and ujab 

habVj^J^KJU- 

A straightforward diagonalization gives the Hamiltonian: 



H2BJJ ^ 2 + Pi + ^2^2 + Pi] , 



(18) 



where, defining A^b = J (tj2 _ ^2^2 ^^2^^ 



2 _ "^+"b-Aab 

— 2 ' "^2 



(19) 
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The operators Qi, Pi and (52, ^2 can be viewed as position and momentum operators of two distinct fictitious particles, 
associated with the modes f and 2, i.e. the Hamiltonian ([TS]) is that of two harmonic oscillators. 

The eigenvalues wi 2 up to order Kf obtained within the Holstein-PrimakofF approach coincide with the zero mode 
frequencies of small amplitude oscillations obtained by the semiclassical approach based on the Gross-Pitaevskii 
equations for the two condensate wave functions in Ref. [l^ (see Eq. (26)) and in Ref. [2^ (see Equation at the 
beginning of Section IV) for the case of equally populated species. When uj\ vanishes, a phase separation takes place, 
resulting in a MQST state. 

Indeed, from Eqs. the stability condition is: 



|Aa.|< J(A, + %) {k, + ^]=KI^, (22) 



where Ajjj^ is the critical value of the interspecies interaction which sets the onset of phase separation regime. Such a 
condition agrees the one given in Ref. [2^ (see Equation (10) in Section IV) and reduces to: 

|Aat| < VAaAb, (23) 

when the limit J" , J** — >■ 00 is taken. 

N „„f , ,2 2 2 , ,2 _ A I S-2 



In the symmetric case Aq = Af, = A, Ka = Kb = K, Na ^ Nt ~ ^ wc get uj^ = uj^ ^ lj^ where uj ~ A^K + K'^, 
and uJab ~ Aah^K. As a consequence A^b = 2aja6 and the eigenvalues (|19p simplify as: 

Wi = - Wafa, Cj| = + Wab , (24) 

which result in the stability condition: 

|A,b|<A + 2^ = A:,. (25) 

In the next Sections we will derive the analytical expressions for the stationary states and discuss the corresponding 
quantum dynamics of the system. 

III. STATIONARY STATES 

Since the Hamiltonian ((T5| is that of two independent particles H = Hi + H2, the corresponding Hilbert space 
is simply given by the tensor product £a ® £b = £1 ® £2 and we can find a basis of eigenvectors for H2BJJ in the 
following form: \ip) = \^^) |</2^), where |(/3^) and |(y5^) are eigenvectors of Hi and H2 within £1 and £2- Since Hi and 
H2 are simply harmonic oscillator Hamiltonians, we could define two pairs of creation and annihilation operators, one 
for each mode, as follows: 



1 

V2 



(26) 



(27) 



being i = 1,2. Now, if we define the ground states of Hi and H2 as \ipQ) and Ifo), we easily obtain eigenvalues and 
eigenvectors within these two subspaces as: 

Ei^{n + ^)huji, |(pi) = -^(a+)'>i), (28) 
E^^^{p+^)h.2, \^l)^^^{atyyo). (29) 
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So the stationary states of the full Hamiltonian (|T8)) are: 

Wn,p) = \ = {alT {<4Y IVo,o) , (30) 



1 \ . f I 



^ \/n\p\ 

and the corresponding energies are: 

En,p ^El + El = + ^ j hui + (^P + f j (31) 

We note that since the Hamiltonian ([T5|) preserved the original parity symmetry of the original Hamiltonian its 
eigenstates could also be classify according to their parity under a — —a and 6 — — &. Since oi and a2 arc linear 
combinations of a, 6, the eigenstates can also be classified by their parity under ai 2 — >■ — ai.2- Using (j30[) . it is then 
clear that the even eigenstates are those with n+p even and the odd eigenstates the ones with n-\- p odd. So we can 
define the parity of a state as (— 1)"+^*. 

We stress that this spectrum is not unbounded because an infinite number of unphysical high energy states have 
been added. Thus a constraint has to be included in order to satisfy the conditions (a^a) ^ 2 J°, {h^h) ^ 2J^ . Solving 
these constraints will give limits to the value of n and p and wc will recover a finite dimensional Hilbert space. Let us 
notice that, through the repeated action of the operators and aj, wc can obtain stationary states of the system 
with a given number of quanta in each mode. The action of a^, ai, aj, a2 on the stationary states \'-Pn,p) is as follows: 

at Wn,p) = \/n+l |(p„+i,p) , ai \'Pn,p) = \/n\Lp„~i,p) (32) 

at Wn,p) = 1 \'fn,i,+ l) , 02 Wn.p) = V^l<^«,p-l) • (33) 

Generically, wi and L02 are incommensurate with each other and there are no degenerate levels since there do not exist 

two different pairs of integers and |7i ,p | such that rtcji + puj2 = n oji + p lu2- Such degeneracy may exist in 

the non-generic case where the ratio ^ is a rational number. In the presence of degeneracy, the non-linear terms that 
we have neglected can lift the degeneracy, unless the states have different parity. 



IV. QUANTUM DYNAMICS 

We are interested in the time evolution of the mean values of the observables J°, J^, that is the population imbalance 
between the left and right well of the potential of each bosonic species. In order to carry out such a program and to 
impose the correct initial conditions it is much more convenient to start from the Heisenberg equations of motion for 
the observables Qi, Q2, Pi, P2- 



d 

di 
d 

di 



(P,) 



1 

ifi 
1 



m,H2BjA) - {P^), 
{[P^,H2BJJ])^~U;f{Q,) 



(34) 
(35) 



which give rise to the following time evolution: 



{P^){t) = (P,) (O)cosa;,i-a;, (g,) (0)sina;,t. 



(36) 

(37) 
(38) 



All we need now is to express J^, in terms of Qi, Q2, Pi, P2 by means of Eqs. (fT4|) . (fT5|) . (|20|) . ([21]) : in this way 
the initial conditions (j^) (0), (j^) (0), (J^) (0), ( J^) (0) are well known. 



Starting from Eqs. 



- (PT|) we find: 



Qi = 



Q2 



J. 



-.X 



b" 



ja 



ja 



(39) 
(40) 



7 



whose inverse transformation gives J" and in terms of Qi,2 and permits us to readily obtain the time-evolution of 
their averages: 



(4) it) 



a 


[Q2){t)-a (Qi) 






a" h' a'b" 











b'^^{Q2) {t)-b" (Qi)(t) 



a o 



a b 



(41) 



(42) 



The coefficients a ,b ,a ,b are defined in the Appendix. The initial conditions relevant for our study are the one 
with a small imbalance between the two wells for each species and the coherent initial states. For the first case we 
choose (J^) (0) = ±1, (4) (0) = ±1, (J^) (0) = 0, (4) (0) = 0, while the particle number is equal to ja = jb = 1000. 
Concerning the chosen values of the interaction strengths Aa, A;, and A^;,, in the following we refer to the mixture of 
^^Rb and ^'^ Rb atoms realized by the JILA group [l6j . 

Figs. [T] and [2] show the dynamics of (J"''') in the case in which there is a small imbalance between the two wells, 
specifically we consider the case in which there is one unit difference in the left and in the right well, in the absence of 
imbalance between the two species (the corresponding parameters are reported in the figure caption). Here we note 
a coherent tunneling between the two wells. 

Figs. [3] and H] show instead the behavior of (J"''') in the case of imbalance between the two species, with an 
imbalance between the two wells of one and two units and for two different values of Aab (0.8 and 1.). As one can 
note, at increasing Aab one approaches a phase separation instability in which the two species tend to separate in the 
different wells. This behavior can be understood in terms of the behavior of the eigenfrequencies uii 2 vs A^;,. In Fig. 
[5] and Fig. [S]one of the two frequency becomes imaginary for a critical value of Aab, thus signalling an instability. 
Let us note that the instability point is a function of A^, Af, and usually takes place for a critical positive value of the 
interspecies interaction, as discussed in Section 2, Eqs. and ([?5|) . In case in which this interaction is attractive 
the system is always in a coherent tunneling regime. 



1.0 
0.5 
0.0 
-0.5 



-l.Ot 



0.0 0.2 0.4 0.6 0.8 



1.0 



t 



FIG. 1: Behavior of the average value of J"'' for Aa = A;, (units of energy), Aab = 2.13, Ka = Kt = 10. and initial conditions 
{JS}{0) = {Jt}{0) = 1-, {JS'''){0) = 0. The time is expressed in units of energy /h. 

A few comments on the dynamics of the system are in order here. Compared to our previous analysis[25|, the 
present analysis does not allow the study of long-time scale phenomena since their detection is abruptly increased 
with N, thus only short-time scale effects are reliable. Furthermore we point out that the dynamics should become 
aperiodic in the general case. 

When the initial state is a coherent spin state for each species, IV' (0)) = IV' (0))a I''/' (*-*))&' where |?/'(0))- = 

E™:l^./2 y^i=fE^tan- (I) e---^- |m.), a = sm^-/^ (f ) cos^'/^ (|) e^^^-, ^ = a,b, then the 
initial conditions are: 

(J«) (0) = -^cosOa, (4) (0) = cos9b, ( Jy") (0) = ^ sin^a sm(j)a, (4) (0) = ^ sui0b sm(f>b. 

In this case the same type of behavior, as for the small imbalance, is observed. In Fig. [7] we take the values 
Oa^Ob = tt/2. and (j)a = tpb = 7r/4. 

The quantum dynamics above investigated could be experimentally reproduced. If we refer for instance to the 
mixture of ^^Rb and ^"^ Rb atoms realized by the JILA group [l^, a wide tuning of s-wave interactions is possible 
via Feshbach resonances. In particular it is possible to fix the scattering length of and to tune the scattering 



8 




t 



FIG. 2: Behavior of the average value of J"'** for Aa = A;, (units of energy), Aab ~ 2.13, Ka = 10., A';, = 10. and initial 
conditions (J°)(0) — 1, {j!^){0) = —1., (Jy'*)(0) = 0. The time is expressed in units of energy/?!. 



2fv 




-1 



-21 1 < ".^ 

0.0 0.1 0.2 0.3 0.4 0.5 

t 

FIG. 3: Behavior of the average value of J"'*" for Aa — At, (units of energy), Aab = 0.8, Ka = 10., Kt = 10. and initial 
conditions (J")(0) = l.,(J^)(0) = 2., {Jy'''){0) = 0. The time is expressed in units of energy/fi. 




t 



FIG. 4: Behavior of the average value of J"'** for Aa = A;, (units of energy), Aa6 = 1., Ka = 10., Kb — 10. and initial conditions 
{Jx){0) = l.,(Ji:)(0) = 2., (Ja'*)(0) = 0. The time is expressed in units of energy/ft. 




FIG. 5: Behavior of tJi(t) (dashed line) and uj2{t) (straight line) for Aa — Ab (units of energy). 
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FIG. 6: Behavior of uji{t) (dashed hne) and ijJ2{t) (straight hne) for Aa — 2Ai, (units of energy). 




t 



FIG. 7: Behavior of the average value of J"'*" for Aa = At (units of energy), A^b — 2.13 and initial conditions with Na — Ni, = 
100. and 6a = Ot = n/2. and <j)a = (pt ~ 7r/4.. The time is expressed in units of energy/?*. 



length of ^^Rb as well as the interspecies one. That allows one to explore the parameter space in a wide range and to 
realize the symmetric regime = A;, = A as well as the asymmetric one. Furthermore one can tune the inter well 
coupling, i. e. the parameters Ka,Kb, in such a way to get the semiclassical limit. Another possible realization of 
the phenomena above described could be obtained with the mixture of '^^K and ^"^Rb atoms produced by the LENS 
group [13], which offers a wide possibility of driving from the weak to the strong interacting regime because of the 
presence of several magnetic Feshbach resonances [33| . 



V. SEMICLASSICAL DYNAMICS 



In this Section we briefly introduce the semiclassical limit of our model within the linear approximation in order 
to make a comparison with the quantum results obtained above. A detailed semiclassical analysis has been already 
carried out in the recent literature (see Refs. [Tsl - E^ ). Here we only recall the classical equations of motion to give a 
physical interpretation of q^.b and Pa^b in Eq. (|14p. From the Hamiltonian we can derive the following equations 

of motion for the components of the vectors: 



-^a,b 

vectors: J : 






rl 


A f 
^^ab^ y ^ z ' 


(43) 


d T'> 

""Jx A jb jb 


A 7° 
^^ab'Jy'J z I 


(44) 


dt ~ " ^ ^ ^ 


AabJ X ^ z ^o.'-^z ' 


(45) 


^ - A J\J'' + 

dt ~ ^ ^ ^ 


AabJ^Jz + ^bJzJ 


(46) 


d 7" 
dt ~ 




(47) 
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J T'' 

-dt ^ ^^^^ 

These equations imply that (T'")^ = (J°)2 + (J^)2 + (ja)2 ^^^^ (j'''-^2 ^ (jh)2 ^ (^6)2 _^ ^jb)2 ^^^.g constants, so we 
can introduce: 

J," = ||J'^||sin^aCOS(^a, Jy" = 1 1 J^l siu 0, siu , = 1 1 J"^ 1 1 COS , (49) 

and: 

4 = ||J''||sin6'6Cosv?f„ = \\j''\\smebSmLpb,J^^ = || cos6'6. (50) 
Using (gni) and ((SDl) in (gS])-®, wc obtain the equations [131 : 

~dt 

— Kbsmipb, (52) 

= ( Aa J" cos +Aa,,j'' cos e^fc )+/;:„ cot 6I„ cos (pa, (53) 

at 

^ = (A6j''cOSeb + AatJ"cOs6'a) +i^bCOt6'tCOS.^h. (54) 

These equations coincide with Eqs. (5)-(8) in Ref. and Eqs. (5) in Ref.[2i| and Eqs. (3) in Ref.[2l|. The 
energy conservation introduces one extra constraint, so that the phase space is actually three-dimensional. This may 
permit in certain conditions the observation of classical chaos. If we linearize the Equations ([?T|) - ([M1) around the 
point Oa = Ob = 7r/2, ^Pa = = 0, we find the equations of motion: 



KaSinipa, (51) 



dt 
dS9b 

dt 
d(Pa 

dt 
dipb 
~dt 



KaVa, (55) 

Kbifib, (56) 

-~iAaJ''S0a + AabJ^SOb) - KJOa, (57) 

-{kbJ''5eb + KhJ^SOa) - KbS9b, (58) 



where 6*0 = 7r/2 + S9a and 9b ^ tt/2 + 59b- These equations derive from the Hamiltonian: 
2 2-1 

H,ff = Kar^ + KbJ"^ + - [{K{rf + Ka.n{59af + {hb{j'f + A^fc j') ((J^a )' + ihab^ j" 59 J9b\ , (59) 

with the Poisson brackets, J''(56'a} = 1 and j''(56'b} = 1. By rcscaling the variable and 59i {i = a,b) as 
ip.^ — > ^jij^ 0i and 59i ^/ J'^Ka59i, we do obtain the corresponding classical hamiltonian of ([T5| . with Poisson 

brackets {ipi, J''59i} = 1. This Hamiltonian can be diagonalized in a standard way by introducing a linear combination 
of the variables (pi and 59i that preserves the Poisson brackets. The diagonalized Hamiltonian will be that of two 
independent classical harmonic oscillators of variables ipi,'P2 and 59i,562. Applying then the Bohr-Sommerfeld 
quantization we do reobtain the spectrum pip , giving the desired connection between the semiclassical and the 
quantum approach. This leads also to a physical interpretation of the conjugate variables ga,fc and Pa,b in Eq- (|14p 

as the azimuthal angles of the pseudospins Tf ' . The full classical solution of Eqs. can be found in Refs. 

MM. 



VI. CONCLUSIONS AND PERSPECTIVES 



In this paper we investigated the quantum dynamics of a Bosc Joscphson junction made of a binary mixture of BECs 
loaded in a double well potential within the two-mode approximation. Focussing on the regime where the number 
of atoms is very large, a mapping onto a SU (2) spin problem together with a Holstein-Primakoff transformation has 



11 



been performed to calculate the time evolution of the imbalance between the two wells. This approach allows one 
to exactly solve the system under the assumption of weak interatomic interactions. The results show an instability 
towards a phase-separation above a critical positive value of the interspecies interaction while the system evolves 
towards a coherent tunneling regime for negative interspecies interactions. The detection of a phase separation could 
be experimentally achieved in current experiments with a mixture of ^^Rb and ^"^Rb atoms (l6j. 

We point out that all the above results are obtained within the linear approximation. It would be interesting to 
extend our model beyond the linear regime; in such a case the classical dynamics may exhibit a chaotic behavior in 
some parameter range because the phase-space is three dimensional. At the quantum level, these features will show up 
in the spectrum as well as the eigenstates of the Hamiltonian. Indeed the Hamiltonian is not time reversal invariant 
because of the terms linear in Jx/Jz, and we expect that the distribution of spacings between energy levels should 
follow the GUE (Gaussian Unitary ensemble) statistics [U- Regarding the dynamics, we conjecture that the short 
time scale behavior of the quantum system will look chaotic, but the long time behavior will not. Such an analysis 
will be carried out in detail in a forthcoming publication. 



Appendix A: Coefficients 



The coefficients a , b , a , b are defined as follows: 



b = 



Uuji 



(Al) 



(A2) 



(A3) 



(A4) 
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